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ABSTRACT 

We have recently proposed a new simplified scenario where blazars are classified 
as flat-spectrum radio quasars (FSRQs) or BL Lacs according to the prescriptions 
of unified schemes, and to a varying combination of Doppler boosted radiation 
from the jet, emission from the accretion disk, the broad line region, and light 
from the host galaxy. Here we extend our approach, previously applied to radio 
and X-ray surveys, to the 7-ray band and, through detailed Monte Carlo sim- 
ulations, compare our predictions to Fermi-LP^T survey data. Our simulations 
are in remarkable agreement with the overall observational results, including the 
percentages of BL Lacs and FSRQs, the fraction of redshift-less objects, and the 
redshift, synchrotron peak, and 7-ray spectral index distributions. The strength 
and large scatter of the oft-debated observed 7-ray - radio flux density correla- 
tion is also reproduced. In addition, we predict that almost 3/4 of Fermi-LAT 
BL Lacs, and basically all of those without redshift determination, are actually 
FSRQs with their emission lines swamped by the non-thermal continuum and 
as such should be considered. Finally, several of the currently unassociated high 
Galactic latitude Fermi sources are expected to be radio-faint blazars displaying 
a pure elliptical galaxy optical spectrum. 

Key words: BL Lacertae objects: general — quasars: emission lines — radiation 
mechanisms: non-thermal — radio continuum: galaxies — gamma-rays: galaxies 
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1 INTRODUCTION 

Blazars are a small subclass of Active Galactic Nuclei 
(AGN) characterized by distinctive and extreme obser- 
vational properties, such as large amplitude and rapid 
variability, superluminal motion, and strong emission 
over the entire electromagnetic spectrum. These peculiar 
sourc es are known since the discovery of QSOs (|Schmidtl 
Il963l '). as 3C273, the first quasar to be discovered, is in 
fact a well known blazar. A major distinction between 
QSOs and blazars is that while the former type of AGN 
mostly emit radiation through accretion onto a supermas- 
sive black hole, blazars also host a jet that is oriented 
closely along the line of sight, within which relativistic 

J iarticles radiate losing their energy in a magnet ic field 
Blandford fc Reed Il978l : lUrrv fc Padovanil Il995l ). This 
jet component is responsible for the non-thermal emis- 
sion from radio to high-energy 7-rays and for the extreme 
properties listed above when the viewing angle is small. 
When this angle is larger than 15 — 20° the relativistic 
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amplification and superluminal effects are much less im- 
portant and jet emission across the electromagnetic spec- 
trum, although still present, may not be dominating the 
flux at IR, optical and X-ray frequencies where other nu- 
clear and non-nuclear components emerge. Under these 
conditions the source is usually called a radio galaxy. 

Although blazars represent a small minority 
AGN, the interest in these sources is growing 
they are being found in increasingly large numbers 
high Galactic latitude surveys perfo r med a t microwave 
and 7-ray energies (fciommi et al.l 20091 : lAbdo et al.l 
l2010al : [Planck Collaboration et aL 2011a ). Blazars repre- 
sent also the most abundant population of extragalactic 
sources at TeV energie^. 

Due to the blend of accretion disk, jet and host 
galaxy emission that is often present in their optical spec- 
trum, the classification of blazars based on observational 
properties is a complex issue. This has resulted in some 
confusion in the literature. 
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In a recent paper lGiommi et al.l l|2012al . hereafter Pa- 
per I) examined the various argumentations and the ex- 
perimental evidence behind the classification of blazars 
that have appeared in the literature, and proposed a 
new scenario, which is based on optical/UV light dilu- 
tion, minimal assumptions on the physical properties of 
the non-thermal-jet emission, and unified schemes. These 
posit that BL Lacs and flat-spectrum radio quasars (FS- 
RQs) are mostly low-excitation (LERGs)/Fanaroff-Riley 
(PR) I and high-excitation (HERGs) /PR II radio galaxies 
with their jets forming a small angle with respect to the 
line of sight (sec Buttigliono ct al. 2010, and references 
therein, for a discussion on LERGs and HERGs). We call 
this new approach the simpUfied blazar view. By means 
of detailed Monte Carlo simulations. Paper I showed that 
this scenario is consistent with the complex observational 
properties of blazars as we know them from all the sur- 
veys carried out so far in the radio and X-ray bands, 
solving at the same time a number of long-standing is- 
sues. 

In this paper we extend the Monte Carlo simulations 
presented in Paper I to the GeV 7-ray band and compare 
our expectations wit h the AGN data in the Fermi-LAT 2- 
yr sou rc e catalogue (iNolan et al.ll20l"^ : lAckermann et al] 
l2011bl : IShaw et al.1 I201I hereafter 2LAC) and the 
Fermi-LAT da,ta of a sample of radio selected bl azars 
l|Giommi et al.ll2012bl : [Planck Collaboration! l2011bl ). We 
note that simulations in the 7-ray band are more com- 
plex than those in the X-ray band since the 7-ray fiux of 
blazars is not due to simple synchrotron or synchrotron- 
self-Compton emission, but is usually attributed to in- 
verse Compton scattering of synchrotron and of external 
photons, or to the superposition of more than one com- 
ponent. Furthermore, the 7-ray flux is a much broader- 
band measurement and is subject to the highest source 
variability. 

Table [T] summarizes the simplified blazar view pro- 
posed in Paper I. Sources with their jets forming small 
angles w.r.t. the line of sight {9 < ^blazar ~ 15 — 20°; 
lUrrv fc Padovanilll995l l. and therefore dominated by non- 
thermal emission, are characterized by low values of the 
equivalent width (EW) and/or the so-called Ca H&K 
break, the latter being a stellar absorption feature found 
in the spectra of elliptical galaxies. However, only LERGs 
belonging to this class are real BL Lacs, that is have in- 
trinsically weak emission lines. HERGs, which are sup- 
posed to have an accretion disk and therefore to display 
strong emission lines, appear to show weak lines only be- 
cause these are swamped by non-thermal emission and 
are therefore "masquerading" BL Lacs {italics in the ta- 
ble) and are in fact misclassified FSRQs. Sources with 
somewhat weaker jet dilution still oriented at small angles 
show some optical features, like the emission lines typi- 
cal of FSRQs. In this category we find also "masquerad- 
ing" radio galaxies, that is sources, which are "bona fide" 
blazars but have their non-thermal emission in the op- 
tical/UV part of the spectrum swamped by the host 
galaxy. Finally, misdirected jets characterize the true ra- 
dio galajcy population. 

In short, starting from truly different objects like 
LERGs and HERGs, the combination of viewing angles 
and strong/weak optical/UV light dilution from the jet. 




0.5 1 
Log (Compton dominance) 

Figure 1. The distribution of the logarithm of Compton Dom- 
inance in the sample of radio selected blazars of I Giommi et aU 
ll2012bD with ly^^^y. < IQl'* Hz. The red solid histogram shows 
the blazars detected by Fermi-hAT, while the dotted his- 
togram represents upper limits for 7-ray undetected blazars. 
The green line is the function that represents the entire pop- 
ulation of radio selected blazars used for our simulations. 



determines the commonly adopted classification as FS- 
RQs, BL Lacs, or radio galaxies. This purely observa- 
tional approach in a number of cases assigns the same 
class to intrinsically different objects. 

Throughout this paper we use a ACDM cosmology 
with Hp = 70 km s~^ Mpc'S ilm = 0.27 and = 0.73 
(|Komatsu et al.ll201ll ). 



2 SIMULATIONS 

Our goal is to estimate the properties of a 7-ray flux- 
limited blazar sample, building upon the simulations pre- 
sented in Paper I. We start by simulating blazars in the 
radio band applying the same ingredients used in Pa- 
per I and then estimate their 7-ray properties. To do so 
we must take into account a series of constraints that 
are directly determined from Fermi-LAT observations of 
well-defined radio selected samples of blazars. Namely: 

(i) the redshift distributions of radio selected and 7- 
ray selected blazars ar e very similar (see Fig. 13 of 
lAckermann et aP l2011bl ) ; moreover the redshift distri- 
butions of i<'ermi-detected and undetected r adio selected 
blaza rs are also not significantly different l|Lister et al.l 
M); 

(ii) a significant fraction (~ 30 — 40%) of radio se- 
lected blazars (mostly FSRQs) are no t detected even in 
the deepest Fermi-L AT observations (iLister et al.ll20l"l] : 
iGiommi et"al]|2012bh : 

(iii) the ratio between inverse Compton and syn- 
chrotron peak luminosities, usually known as Compton 
Dominance (CD), is crucial to p redict the 7-ray inten - 
sity of a blazar from its radio flux. ICiommi et al.l (|2012lJ ) 
showed that the CD distribution of radio selected blazars 
is very different below and above i^peak — 10^* Hz, where 
'^pcak is the peak frequency of the synchrotron pow er. Fig. 
□ (adapted from Fig. 22 of lCiommi et al.l l|2012bl 'l) shows 
this distribution for blazars with t'peak < 10^* Hz. The 
red solid histogram refers to all radio selected blazars 
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Table 1. The simplified blazar view scenario. See text for details. 



LERG 



HERG 



viewing angle 



strong jet dilution (EW < 5 A; Ca H&K < 0.4) BL Lac BL Lac^'^^ 8 < ^biazar 

weak jet dilution radio galaxy^^^ FSRQ 6 < ^blazar 

misdirected jet radio galaxy radio galaxy 8 > Sblazar 

Italics denote "masquerading" sources (see text), misclassificd FSRQ, misclassified BL Lac 



that were detected by Fermi-LAT. The dotted histogram 
refers instead to upper limits derived from the Spectral 
Energy Distribution (SED) of those blazars that were not 
dete cted even after 27 mon ths of Fermi-LAT observations 
(see[G iommi et al.ll2012bl . for details). These undetected 
7-ray sources represent a significant fraction of the sam- 
ple and certainly cannot be ignored. The green line is our 
educated guess of the underlying intrinsic distribution. 



> 10^'' Hz (which are all 



The CDs of blazars with fpeak 
detected in the 7- ray band) are instead clustered around 
~ 1 (see Fig. 22 of lGiommi et al.ir2012bl ). We assume here 
that their distribution is represented by a gaussian with 
(CD) = 0.9 and a = 0.4. 

(iv) the sensitivity of Fermi-LAT is a strong func- 
tion of the 7-ray spectr al index F (see, e.g.. Fig. 14 of 
lAckermann et aul2011bl ). We therefore need to to esti- 
mate the latter parameter as well. Since we know that 
for b lazars this correlate s strongly wit h u^^^^, (e.g., Fig. 
29 oflAbdo et al.l (12010^ ) and Fig. 17 of lAckermann et al.l 
(|2011bl )). we use that depen dence for our purposes. By 
fitting the data of Fig. 17 of Ac kermann et al.l (|2011bh we 
derive F — —0.22 x log{Up^^^) + 5.25, with a dispersion 
a = 0.15. 



2.1 The dependence of 7-ray flux on radio flux 
density 

The dependence of 7-ray flux on radio flux den- 
sity has been discusse d for quite some time (e.g., 
lAckermann et aDl2011al . and references therein). In this 
paper we will base all our calculations on the observa- 
tional constraints listed above. Based on point (i) above 
we therefore assume that the 7-ray to radio flux density 
ratio, f-y/fr, is independent of redshift, and, in line with 
the minimal assumption approach of Paper I (and con- 
sistently with observational data), we assume that this 
ratio is also independent of luminosity; hence f-y/fr — 
F{CD, 



Since the CD distribution depends on 



for 

'^peak < 10^* Hz we use a linear dependence on CD 
with parameters derived from the sample of radio selected 
blazars: log .f~, / fzoHz = 0.5 x log(CD) — 8.0 (where 
is in units ph cm~^ s~^ and fsoHz is in Jansky). For 
higher I'peak values, where inverse Compton scattering 
occurs in the Klein-Nishina regime and all sources have 
CD ~ 1, we use instead a correlation between the 7- 
ray to radio flux ratio and z-'pcaki that is log / fscHz ~ 
0.13 xlog(;.pS,k)- 9.8. 



2.2 Simulation steps 

Our Monte Carlo simulations consist of the sequential 
execution of the following steps: 



(i) all those described in detail in Paper I to simulate 
the radio through X-ray SED of blazars. In short, these 
include (readers familiar with Paper I can skip directly 
to point (ii)): a) drawing values of the radio luminosity 
and redshift based on our assumed luminosity function 
(LF) and evolution (see below); b) drawing a value of the 
Lorentz factor of the electrons radiating at the peak of the 
synchrotron power and, based on that, calculate fpeak) 
suming a simple synchrotron self-Compton model (with 
B = 0.15 Gauss and a Doppler factor S randomly drawn 
from a gaussian distribution with (5) = 15 and a— 2); 
c) calculating the observed radio flux density and the 
non-thermal emission in the optical and X-ray bands; d) 
adding an accretion component (only for beamed FR II 
sources), re-scaling an SDSS quasar template to a value, 
which depends on radio power; e) drawing a value of the 
EW of Lya, C IV, C HI, Mg II, H/3, Ha starting from the 
EW distribution of the broad lines of SDSS radio quiet 
QSOs; f) adding the optical light of the host galaxy as- 
suming a standard giant elliptical as observed in blazars; 
g) calculating the total optical light and the observed EW 
of all the broad lines by taking into account the dilution 
due to the non-thermal and host galaxy optical light. In 
order to keep our simulations as simple as possible, we 
used no joint probability distributions during the pro- 
cess. In fact, as already discussed in Paper I, some of the 
observed correlation between different blazar quantities 
could be well reproduced as selection effects in a simple 
scenario with no intrinsic correlations; 

(ii) as regards the LF, since the radio powers of 7- 
ray-selected sources reach lower values than the radio- 
selected ones, we extrapolated the radio LF of Paper I, 
that is >I>(P) (X (in units of Gpc~^ P~^) by a factor 
~ 4 down to 5.0 x 10^'' W/Hz assuming the same slope. 
As for the evolution, we use the same model of the type 
P{z) = (H- z)'=+'^^ which allows for a maximum in the 
luminosity evolution followed by a decline, with k — 7.0 
and P = —1.5 (which implies a peak at 2: ~ 1.78). This is 
only slightly different (within la) from the k value used 
in Paper I (fc = 7.3, while /3 is the same) but it gives 
a better agreement with the 7-ray observables without 
affecting the results of Paper I, as discussed in Sect. 5.4 
of that paper. 

As detailed in Paper I, we take into account the well- 
known observational result that low-power radio sources 
display a much weaker cosmological evolution than high- 
power ones by having a fraction of non-evolving sources 
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Gamma-ray flux[E>1 OOMeV]{ph cm-^ s-^) 



Figure 2. Tlie spectral index of blazars in our simulated sur- 
vey of 1,000 7-ray sources is plotted against 7-ray flux. FSRQs 
are plotted as filled red circles, BL Lacs as open blue circles. 
The flux cut applied reproduces the sensitivi ty limits of Fermi- 
LAT i n the 2LAC catalogue (see Fig. 14 of lAckermann et al.l 
l2011hl V 

equal to 1 for Pr ^ 5 x 10^* W Hz~^ and decreasing 
monotonically with power to for Pr ^ 5x 10^^ W Hz^^; 

(iii) for sources with t'peak < 10^^ Hz draw a value of 
the CD from the distribution shown as a green line in Fig. 
[T]and then use the f^/fr - CD correlation, while for those 
with I'peak > 10^* Hz use the f-y/fr - t'poak correlation 
fSect. [2TT|) . Get by multiplying the resulting f-y/fscHz 
by the radio flux density calculated above; 

(iv) obtain a value of the 7-ray spectral index from the 
correlation with i^p^ak described in Sect. [5] 

(v) apply the Fermi 2LAC sensitivity limits, which are 
a strong f unction of the 7-ray spectr al index as shown in 
Fig. 14 of lAckermann et al.l (|2011bl '). The same plot for 
our simulation is shown in Fig. [2] 

As done in Paper I, we classify a source as an FSRQ 
if the rest-frame EW of at least one of the broad lines 
that enter the optical band in the observer frame (which 
we assume to cover the 3, 800 — 8, 000 A range) is > 5 A. 
Otherwise, the object is classified as a BL Lac, unless 
the host galaxy dominates the optical light causing the 
Ca H&K break to be larger than 0.4, in which case the 
source is classified as a radio galaxy. A BL Lac whose 
maximum EW is < 2 A, or for which the non-thermal 
light is at least a factor 10 larger then that of the host 
galaxy, is deemed to have a redshift which cannot be 
typically measured. 



3 RESULTS: COMPARING SIMULATIONS 
WITH Fermi 7-ray DATA 

We now make a detailed comparison between the results 
of our simulations an d the observational data from the 
2LAC "clean sample" (jAckermann et al]|2011bl l in terms 
of fractions of BL Lacs and FSRQs, and redshift, z^p^aki 
and 7-ray spectral index means and distributions. To in- 
clude the latest measurements we have updated the op- 
tical classifications and redshift estimates of the "clean 
sa mple" using t he rec ently published spectroscopy results 
of|Sh aw et al ] l|2013l ). As stressed in Paper I, the aim of 
our simulations is not to reproduce the fine observational 



Table 2. Results from a simulation of a 7-ray survey 




Source type 


Number of 
sources* 




(1°S ''leak)* 


(r>* 


FSRQs 
BL Lacs 
Radio galaxies 


349.4 

549.6 (277.3)" 
101.0 


1.23 
0.58 
0.06 


13.15 
14.99 
12.98 


2.34 
1.89 
2.39 


Total 


1,000.0 


0.82 


14.15 


2.10 


'Average of 10 runs each simulating 1,000 sources 
**BL Lacs with measurable redshift 




Table 3. 2LAC data 








Source type 


Number of 
sources 




(l°g ''leak) 


(r) 


FSRQs 
BL Lacs 
Unclass. blazars 


309 

418 (168)** 
128 (14)** 


1.18 
0.44 
0.32 


13.23 (72%*) 
15.05 (76%*) 
15.14 (52%*) 


2.36 
1.98 
2.07 



* Percentage of sources with measured w 
** Sources with measured redshift 



details. Our approach is instead to reach robust conclu- 
sions keeping the number of assumptions to a minimum. 
Our simulations include 1,000 7-ray detected blazars, a 
number similar to that observed in the 2LAC catalog. To 
ensure statistical stability of the results the simulation 
was repeated ten times and the results were averaged. 

Table [2] summarizes our main results by giving the 
number of sources per class and their mean redshift, 
^'peaki 7-ray spectral index. The number in paren- 
thesis refers to the BL Lacs with measurable redshift, 
to which the mean redshift pertain. Table[3]gives instead 
the observed values for the 2LAC sample, where the num- 
bers in parenthesis in the !^|oak column give the fraction 
of sources for which this parameter could be estimated. 

About 61% of our blazars are classified as BL Lacs. 
This agrees well with the observed value, which ranges be- 
tween 56±4% and 64 ± 4% (where the la errors are based 
on binomial statistics: lGehrelslll986l ). The first number 
is derived by ass u ming t hat the unclassified blazars in 
lAckermann et all l|2011bl ) are all FSRQs while the sec- 
ond one presumes that they are all BL Lacs. Given that 
their mean i^|,ak S'ncl 7-ray spectral index are closer to 
those of BL Lacs, the second value appears more realistic. 

The 2LAC mean value is also in good agreement 
with the mean redshift for our simulated FSRQs (1.18 
vs. 1.23), while for BL Lacs it is somewhat smaller (0.44 
vs. 0.58). Fig. [3] shows the overall agreement between 
our simulated redshift distributions (where we have only 
included sources with a measurable redshift) and the ob- 
served ones. 

70% of the BL Lacs (44% of those with redshift) have 
a standard accretion disk and are therefore FSRQs, which 
are classified as BL Lacs only because their observable 
emission lines are swamped by the non-thermal contin- 
uum, as was the case for many radio-selected blazars in 
Paper I. 

About 50% of our BL Lacs have a redshift deter- 
min ation, which is cl ose to the 2LA C value of 40 ± 
4% l|Shawet all 120131 '). As stated in lAckermann et aP 
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Figure 3. The redshift distribution of t he blazars in the 
2LAC catalogue llAckermann et alJlioTTbl : IShaw et al.ll2013h 
(black solid histogram) and that of our simulations (red dot- 
dashed histogram) for FSRQs (top panel) and BL Lacs (bot- 
tom panel). 




y-ray photon spectral index 



Figure 5. The 7-ray spectral (photon) index distribution of 
FSRQs (top panel) and BL Lacs (bottom panel). Black solid 
and red dot-dashed histograms represent 2LAC and simulated 
blazars respectively. 




Figure 4. The distribution of synchrotron peak energies for 
blazars in the 2LAC catalogue (black solid histogram) and 
that of our simulations (red dot-dashed histogram) for FSRQs 
(top panel) and BL Lacs (bottom panel). 



(|201Hj ). many of the unclassified blazars do not have 
an optical spectrum so we cannot include them in the 
evaluation of the fraction of BL Lacs for which a redshift 
could not be obtained because they were totally feature- 
less. According to our simulations, the "real" redshift of 
featureless BL Lacs ranges typically between 0.5 and 2.5, 
which is much larger than that of those with redshift. Ba- 
sically all of these sources possess an accretion disk and 
therefore are FSRQs. 

A small fraction (~ 10%) of the simulated blazars 
are classified as radio galaxies, again as was the case in 
Paper I. These are bona-fide blazars misclassified because 
their non-thermal radiation is not strong enough to dilute 
the host galaxy component. Their relevance is discussed 
in Sect. 1431 

The predictions for t'pcak agree well with observa- 
tions, as shown in Tab.[5]and[3] The 2LAC mean value is 
(log J^peak) = 15.07 if one includes the unclassified blazars 
with the BL Lacs. Fig. |4] compares the distributions of 



'^pcak for sources classified as FSRQs and BL Lacs in our 
simulation with those of the 2LAC blazars. The agree- 
ment is clearly quite good and reproduces the well known 
fact that 7-ray-selected BL Lacs ten d to have val- 
ues significantly high er than FSRQs (|Abdo et al.ll2010bl : 
iGiommi et al.|[2012bl ). One caveat here is that only about 
3/4 of 2LAC blazars have Vp^zk values, which moreover 
were not estimated from a fit to the SED but by using 
the radio-optical and X-ray-optical broadband indices. 
The 

^pcak values for our simulated BL Lacs without 
redshift peak at ~ 10^* Hz, as for the 2LAC sources. 
Only ~ 10% of simulated BL Lacs are Low Synchrotron 
Peaked (LSP: u^^^^ < lO" Hz), with ~ 33% being In- 
termediate Synchrotron Peaked (ISP: 10^* Hz < /^p^ak < 
lO^'^ Hz) and ~ 57% High Synchrotron Peaked (HSP: 
'^pcak > 10^^ Hz). The corresponding 2LAC values are 
~ 20%, ~ 25%, and ~ 55%, which are not too different, 
considering also the caveat mentioned above. Very simi- 
lar fractions (~ 23%, ~ 20%, and ~ 57%) are obtained 
by including with the BL Lacs the unclassified blazars. 

Finally, the 7-ray spectral photon indexes are also 
consistent with observations, with a 2LAC mean value of 
(r) = 2.00 if one includes the unclassified blazars with 
the BL Lacs. Fig. [5] compares the distributions for our 
simulated FSRQs and BL Lacs with those of the 2LAC 
blazars. Despite the slight shift in the peak of the distri- 
bution for BL Lacs, the agreement is again quite good 
and reproduces the well known fact that 7-ray-selected 
BL Lacs tend to have 7-ray s pectral indexes fi atter than 
those of FSRQs (e.g. lAckermann et al.ll2011bl . and refer- 
ences therein). 

3.1 7-ray predictions for radio selected blazars 

To further check our scenario and verify that we are con- 
sistent with the observational constraints listed in Sect. 
[2] we have simulated a sample of radio selected blazars 
with a fiux density limit of 1 Jy. We have then compared 
our 7-ray predictions with the Fermi-LAT results for the 
sample of radio bright blazars reported in lGiommi et al.l 
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l|2012bl '). We find that about 71% and 91% of simulated 
radio selected FSRQs and BL Lacs respectively are above 
the Fermi 2LAC sensitivity limit. These percentages com- 
pare very well with the observed Fermi-LAT detection 
rate for FSRQs (Jlt\l%) and BL Lacs (100+33%) in the 
radi o sample. 

lAckermann et al] l|2011bl ) showed that the redshift 
distribution of 7-ray selected blazars is very similar to 
that of bright radio selected blazars (see their Fig. 13 for 
a comparison of the 2LAC and the WMAP blazar sam- 
ples). The redshift distributions of our simulated radio 
and 7-ray selected blazars are also very similar with av- 
erage values close to the observed values {z)radio ~ 1-2 
and {z)^-ray = 1.2 for FSRQs, and {z)radio = 0.7 and 
{z)^-ray ~ 0.6 for BL Lacs with redshift determination. 
We make one more comparison between simulations and 
observations in Sect. 14.31 

Given the similarities in the redshift distributions 
and considering that the 7-ray to radio flux density ra- 
tio is independent of redshift, we expect the cosmolog- 
ical evolutionary properties of radio and 7-ray selected 
blazars to be the same. 

3.2 Assessing the stability of our results 

As we did in Paper 1 we assessed the dependence of 
our results on the adopted LF and evolution by mak- 
ing two checks: first, we varied their input values by la; 
second, we used as an alternative LF the sum of the 
B L Lac and FSRQ LFs b ased on the beaming model 
of lUrrv fc Padovanil (|l995l ') (converted to Ho = 70 km 

Mpc~^) assuming, for consistency with the way they 
were derived, a pure luminosity evolution of the type 
P(z) = P{0)exp[T(z)/T], where T{z) is the look-back 
time, and r = 0.33 (consistent with the evolution of 
DXRBS F SRQs and BL Lacs combined, based on the 
samples in lPadovani et aDl2007l '). We also tested the sen- 
sitivity of our simulations to different assumptions on the 
Doppler factors using distributions with mean values that 
ranged between 5 and 30 depending on the radio lumi- 
nosity of the simulated blazar. 

In all cases our main results, that is the prevalence 
of BL Lacs in 7-ray selected samples and the higher red- 
shifts, lower i^pcaki s-nd larger 7-ray spectral indexes of 
FSRQs were confirmed, which shows that they are inde- 
pendent of the details of the LF, evolution, and assumed 
Doppler factor. 



4 DISCUSSION 

4.1 Fermi BL Lacs without redshift 

Despite extensive spectroscopic observations 56 ± 5% of 
the BL Lacs in the 2LAC do not have a redshift de- 
terminatio n due to the lack of feat ures in their optical 
spectrum |Ackermann et al]|2011bh . This high fraction 
of redshift-less blazars is matched very well by our sim- 
ulations, which indeed predict a value ~ 50%. Most of 
these sources are of the HSP type, ^1/3 are ISP, while 
only ~ 1/10 are LSP. These fractions are also in agree- 
ment with what is observed in the 2LAC sample. 




l-og(L„„,_,o„s,v|)[erg s-'] 

Figure 6. The 7-ray spectral index vs. 7-ray luminosity plot 
for 7-ray flux limited samples; FSRQs: filled red circles, BL 
Lacs: open blue circles. Top panel: the case of the 2LAC cat- 
alogue jAckerman n et al. 2011b) ; middle panel: simulated 7- 
ray flux limited survey where objects with no redshift infor- 
mation are not plotted; bottom panel: simulated 7-ray sur- 
vey where the blazars with no redshift determination appear 
as gray open squares and sources classified as radio galaxies 
appear as light green symbols. The dashed line gives the ap- 
proximate boundary between LERG and HERG blazars in the 
simulation. 



Our simulations imply that these sources have much 
larger redshifts than those of BL Lacs with redshift in- 
formation: while in the latter case ~ 62% of objects have 
redshift ^0.5 , in the former case redshift values range 
mostly between 0.5 and 2.5, with peak probability be- 
tween 0.9 and 1.2. Therefore, one canno t assume for them 
a valu e typical of BL Lacs (as done bv lAckermann et al.l 
bOllbl ') but larger values (2 ~ 1) need to be used. More- 
over, most of these sources are predicted to be quasars 
with their emission lines heavily diluted by the non- 
thermal continuum. These should therefore be included 
with the FSRQs when studying number counts, cosmo- 
logical evolution, and LFs, since their exclusion will bias 
the results. 

We note that the existence of such sources has 

been recently confirmed at relatively high redshift. 

IPadovani. Giommi. fc Raul (l2012l ) have studied the SEDs 
of eleven 2LAC blazars without spectroscopic redshift 
but for which a photometric redshift > 1.2 was derived 
by fitting SED templates to their UV-to-near-IR multi- 
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Figure 7. The 7-ray spectral index vs. 7-ray luminosity 
plot for radio flux limited samples; FSRQs: filled red cir- 
cles, BL Lacs: open blue circles. Top panel: the case of the 
Planck Swift Fermi radio b right sample (data taken from 
IPlanck Collaboration! 12011^ 1: middle panel: simulated radio 
flux limited survey where objects with 7-ray flux below the 
2LAC catalogue limit and those with no redshift information 
are not plotted; bottom panel: simulated radio survey where 
the blazars with 7-ray flux below the 2LAC limits appear as 
light green points and sources with no redshift determination 
appear as gray open squares. 



band photometry obtained q uasi-simultaneousl y with 
Swift/\]VOT and GROND by iRau et all (|2012D . Four 
blazars turned out to be of a type never seen before but 
with properties expected by us for some of the BL Lacs 
without spectroscopic redshift: large i^pcak (~ 5 x 10^^ 
Hz) and high-power (Lpoak « 6x lO**® erg s~^). Four more 
were ISP and two were LSP, all six with Lpcak ^ 2 x 10*'' 
erg s"'^, which is well into the FSRQ regime. It then looks 
like these sources are very plausible FSRQs in disguise. 



4.2 "Masquerading" BL Lacs 

How can one, in practice and in general, spot the "mas- 
querading" BL Lacs? There are a few options: 

(i) if a photometric redshift esti mate is ava ilable, as in 
the case of the sources studied bvlRau et all (2012), one 
can derive the radio, or synchrotron peak, power and see 
if they are FSRQ-Uke; 

(ii) if infrared spectroscopy is available then there is 
a good chance that the source could turn out to be an 



FSRQ. Our simulations of the 2LAC, in fact, imply that 
~ 50% of the BL Lacs without redshift should show a 
strong (EW > 5 A) Ha line, which for 2: ~ 1 would fall, 
for example, at ~ 1.3 yLxn. For a further ~ 25% one should 
find 2 < EW < 5 A and therefore it should be possible 
to assign a redshift to these sources as well; 

(iii) finally, as mentioned in Sect. [31 basically all 
sources without redshift in the 2LAC catalogue are pre- 
dicted by our simulations to have an accretion disk and 
therefore to be intrinsically quasars. So for this sample 
the mere fact of having a featureless spectrum should be 
enough for a source to be considered, perhaps ironically, 
an FSRQ with high probability. We stress that this result 
is sample- and fiux density limit-dependent. 

In summary, we notice that in our scenario the his- 
torical 5 A separation between BL Lacs and FSRQs has 
no physical meaning and that the detection of any line as- 
sociated with the broad line region of any strength should 
be sufficient to grant a source FSRQ status. 



4.3 The 7-ray - spectral slope luminosity plane 

iGhisellini et"al] (|201ll , |2012| ) noted a strong correlation 
between the lOOMeV - lOGeV luminosity and the 7-ray 
spectral slope in the 2LAC sample, arguing that this re- 
lationship reflects a physical distinction between FSRQs 
and BL Lacs. In this interpretation the presence (or ab- 
sence) of strong emission lines is associated to a regime of 
high (low) radiating efficiency and, consequently, larger 
(or smaller) 7-ray luminosity. 

We tested if such a correlation could be induced 
by the Fermi-LAT sensitivity dependence on spectral 
slope by comparing the 2LAC data with our simulated 
7-ray survey blazars. Fig. [6] plots the 7-ray slope vs. 
lOOMeV - lOGeV luminosity for blazars with a mea- 
surable redshift in the 2LAC clean sample (top panel), 
and in our simulated survey (central panel). The trend 
of high-luminosity /steep 7-ray spectrum associated to 
FSRQs (filled red points) and low luminosity /flat 7- 
ray spectrum associated to BL Lacs (open blue points) 
is clearly present also in the simulation, with BL Lacs 
mostly concentrated in the flat-spectrum/low luminos- 
ity corner and FSRQs populating almost exclusively the 
steep spectrum/high luminosity part of the plot. If we 
consider also BL Lacs with no redshift determination and 
objects (mis)classifled as "radio galaxies" (gray squares 
and fight green open circles respectively in the bottom 
panel of Fig. (G)) we see that the 7-ray - spectral slope lu- 
minosity plane is filled in all its quadrants, although the 
different densities imply that a strong correlation is still 
present. We note that below ~ 10** erg s~^ we predict 
the existence of BL Lacs with photon indexes somewhat 
steeper than those observed in the 2LAC sample. This 
is likely due to our simple power law extrapolation of 
the radio LF and to the misclassification of some sources 
close to the BL Lac/ "radio galaxy" border. The LF ex- 
trapolation is also behind the somewhat higher fraction 
of simulated sources with ij 10*^ erg s^^ as compared to 
the observed one. 

We then c onclu d e tha t the trend reported by 
IGhisellini et~aD (|201ll . I2OI2I ) comes out naturally from 
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our simulations and reflects the complex Fermi-LAT sen- 
sitivity limits. However, we also note that the separation 
between "real" (that is, LERG) BL Lacs and FSRQs, 
that in our simulation fall to the left and to the right of 
the das hed line of Fig. [6l is in a greement with the deduc- 
tions of iGhisellini et ail l|2012l ) for sources with known 
redshift. In our scenario BL Lacs with no redshift are al- 
most all expected to be to the right of the dashed line, 
and therefore to be HERGs and diluted FSRQs. 

To better understand the 7-ray spectral slope - lu- 
minosity plane we simulated the 7-ray properties of a 
sample of blazars with a radio flux density limit equal to 
those o f the sample considered bv iPlanck CollaboratiorJ 
l|2011bl ) and compared our exp ectations with t he Fer mi- 
LAT measurement reported bv lGiommi et all (|2012bl ). 

Fig. [7] is the same as Fig. [6] for the case of a radio 
flux density limited sample (see Sect. 13.1]) . The top panel 
plots the 7-ray spect ral slope versus luminosity fo r the 
sample described in [ Planck Collaborationl (|2011ljl and 
I Giommi et al.l (|2012bl ). The left pointing arrows repre- 
sent blazars not detected even in the 27 month Fermi- 
LAT data to which a spectral slope of 2.5 was assigned. 
BL Lacs with no redshift determination are not plotted. 
The results of our simulations are shown in the middle 
and lower panels where we plot respectively FSRQs and 
BL Lacs expected to be detected by _Fermi-LAT in a 27 
month integration as red and blue points, while BL Lacs 
with no redshift determination are shown as gray open 
symbols, and blazars undetected by _Fermi-LAT appear 
as light green circles. 

As apparent from Fig.[7]in this radio selected sample 
of blazars the correlatio n between 7-ray spectral slo pe 
and luminosity noted by IGhisellini et al.l (|201ll . l2012h is 
not present, in good agreement with our simulations. This 
confirms that it is crucial to fully understand the biases 
of the various surveys before one can draw any conclusion 
on the underlying physics. 



4.4 The blazar 7-ray - radio flux density 
correlation 

Correlations between the blazar 7-ray and radio fiux 
densities have been discussed by several authors in the 
past few years without reaching definite conclusions. This 
was due to limitations connected to the availability of 
large and unbiased samples, lack of simultaneity, and 
presence of a significant number of non-detections. Re- 
cently, a consensus seems to have been reached on the 
existence of a correlation bet ween the two bands, al- 



beit with a large scatter (e.g., Ackermann et af] l2011al : 



I Giommi et al.l l2012bl : iLeon-Tavares'et'^ 20121 . and ref- 
erences therein). 

Fig. [8] shows the 7-ray fiux vs radio flux density for 
our simulated 7-ray survey. This is remarkably similar to 
the equivalent plot base d on Fermt-hPC£ data (Fig. 3 of 
lAckermann et al]|2011al ). 

Despite the large spread in f-f/fr, induced by the 
wide dynamical range of the observed CD distribution 
(which can result in 7-ray fluxes up to a factor 100 or 
more different for the same radio flux density) , a correla- 
tion is clearly present in the 2LAC simulation shown in 




1000 

Radio flux density [mjy] 

Figure 8. Tiie 7-ray flux vs. radio flux density in our 2LAC 
simulated survey (reduced to 1,000 sources for better visibil- 
ity). Fifled red circles represent FSRQs, open blue circles rep- 
resent BL Lacs. The dashed line corresponds to the 7-ray flux 
limit of a survey a factor fs! 5 shallower than 2LAC. 



Fig. [8l For 7-ray surveys with much shorter integration 
times (e.g. corresponding to the dashed line in Fig. [8|, 
the 7-ray dynamical range decreases significantly while 
the radio one remains similar. This, combined with the 
much reduced statistics, makes it much more difficult to 
detect the intrinsic 7-ray-radio correlation. 



4.5 Fermi-LAT unassociated high-Galactic 
latitude sources 

As mentioned in Sect. O the scenario we propose pre- 
dicts the existence of "masquerading" radio galaxies, 
that is blazars having their optical non-thermal emission 
swamped by the host galaxy. These are intrinsically rare 
sources with high CD values that, when detected close to 
the fermi-LAT sensitivity limit, have a radio and opti- 
cal non-thermal component significantly lower than that 
of blazars with smaller CD. They would then display a 
pure elliptical galaxy optical spectrum and would appear 
as radio sources typically fainter than currently known 
Fermi steep spectrum blazars (that are also character- 
ized by large LAT error regi ons) thereby failin g the as - 
sociation criteria applied bv lAckermann et al.l l|2011bl ). 
Objects of this type could then be the counterparts of a 
significant fraction of the still unassociated high-latitude 
Fermi 7-ray sources, which make up ~ 20% of the total 
at |fo| > 10°. This is the same mechanism at work in shal- 
low X-ray surveys where intrinsically rare objects with 
high fx/ fr, that is HSP, are preferentially selected over 
the more common LSP. 

A systematic search for such elusive counterparts of 
Fermi-\jPC£ sources is beyond the scope of this paper and 
will be carried out in the future. Fig. [9] shows, as an 
illustrative example, the SED of SDSSJ173130. 8-^5429, 
which might be a source of this type within the error 
ellipse of the presently unassociated 7-ray source 2FGL 
J1730.8-f5427. 
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Figure 9. The SED of the radio source coincident with the 
17 magnitude z= 0.238 elliptical galaxy SDSSJ173130.8+5429 
within the error ellipse of the unassociated Fermi 7-ray source 
2FGL J1730.8+5427. An example of a high CD radio galaxy as 
a possible counterpart of unassociated _Fermi-LAT high galac- 
tic latitude sources. The green line represents the template of 
a giant elliptical galaxy llMannucci et al]|200lh . 



5 CONCLUSIONS 

We have carried out Monte Carlo simulations of the ra- 
dio through 7-ray emission of blazars selected in the 7- 
ray and radio bands within the framework of the newly 
proposed simplified blazar view presented in Paper I. We 
have extended these simulations to the 7-ray domain us- 
ing a 7-ray to radio flux density ratio that is constrained 
by the relation s hips be tween SED parameters discovered 
bv lAbdo et alJ (|2010bl ') an d by the distribution o f Comp- 
ton Dominance derived bv lGiommi et all (|2012bl ). No de- 
pendence on redshift or luminosity was assumed. 
Our main results can be summarized as follows: 

(1) Our results match very well the blazar content and 
the observational prope rties of the 7-ray sources included 
in the 2LAC catalogue (|Ackermann et al]|2011lj ). In par- 
ticular, we reproduce all the main features of the 2LAC 
sample, namely the prevalence of BL Lacs, the larger red- 
shifts and 7-ray spectral indices and lower i^poak of FS- 
RQs as compared to BL Lacs, and the fraction of BL Lacs 
without redshift, including their split in terms of fpe ak • 

(2) The simulations also reproduce well the 7-ray 
propert ies of the radio flux density lim ited sample consid- 
ered by lPlanck Collaboration! l|2011bl ) and|G iommi et al.l 
(|2012bh . 

(3) Our simulations imply that the redshift of the BL 
Lacs with featureless optical spectrum is much larger 
than that of those with measured redshift, with values 
ranging from 0.5 to 2 and peak probability around 1. 
These sources are also predicted to possess an accretion 
disk and therefore not to be real BL Lacs but quasars, 
and as such should be considered. 

(4) The tendency for BL Lacs of appearing as low 
7-ray luminosity /flat spectrum sources and of FSRQs 
of being hi gh luminosity /steep 7-ray sources noted by 
iGhisellinT et aD ((201 ll I2OI2I ) is also present m our smiu- 
lated 2LAC survey, but not in real and simulated radio 
flux-density limited samples. In both our scenari o and in 
the interpretation of lChisellini et al.l l|201ll . I2OI2I I low lu- 



minosity BL Lacs are low accretion LERG objects, while 
FSRQs and high luminosity steep spectrum BL Lacs are 
HERG sources; ho wever in our simulations, con trary to 
the assumptions of lChisellini et al] l|201ll . [2012l l. almost 
all the redshift-less BL Lacs are diluted FSRQs. 

(5) The strength and large scatter of the observed 7- 
ray - radio flux density correlation, which has caused 
some discussion in the literature, are reproduced by our 
simulations. The latter is induced by the wide dynamical 
range of the observed Compton Dominance distribution. 

(6) About 10% of the 7-ray sources are classified as 
radio-galaxies on the basis of their optical properties but 
are instead moderately beamed blazars with their non- 
thermal emission swamped by the host galaxy. These ob- 
jects fail to pass the association criteria typically applied 
to 7-ray samples and could therefore explain a signifi- 
cant fraction of the still unassociated high-latitude Fermi 
sources. 
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